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Abstract: 

Deep inelastic scattering (DIS) events, selected from 1993 data taken by the 
HI experiment at HERA, are studied in the Breit frame of reference. The frag- 
mentation function of the quark is compared with those of e + e~ data. It is shown 
that certain aspects of the quarks emerging from within the proton in e~p interac- 
tions are essentially the same as those of quarks pair-created from the vacuum in 
e + e~ annihilation. The measured area, peak position and width of the fragmenta- 
tion function show that the kinematic evolution variable, equivalent to the e + e~ 
squared centre of mass energy, is in the Breit frame the invariant square of the 
four-momentum transfer. We comment on the extent to which we have evidence 
for coherence effects in parton showers. 
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1. Introduction 

This paper analyses the spectra of hadrons produced in deep inelastic scattering (DIS) e~p 
collisions as a function of both the four-momentum squared (q 2 = —Q 2 ) transferred from the 
electron and the Bjorken scaling variable, x. In terms of the proton four momentum P, x 
is defined as Q 2 /2P ■ q, but it may be naively thought of as the momentum fraction of the 
struck quark in the proton. We also select on the basis of the dimensionless inelasticity variable 
y = Q 2 / xs where i/i is the overall e~p centre of mass energy. The HI collaboration has already 
published charged particle spectra, analysed in the hadronic centre of mass (CMS) frame, in 
intervals of the total hadronic mass, W (where W 2 ~ Q 2 (l/x — 1)), using 22.5 nb _1 of data 
obtained in the first year of HERA operation U . It is the aim of the present analysis to use the 
increased statistics of 1993 data to cast such spectra into the form of fragmentation functions 
which are reliably related to the hadronisation of quarks. Aspects of these distributions can 
then be used in a comparison with corresponding data from e + e~ experiments |1 to determine 
the relevant evolution variable for such distributions. 

Since the initial state of e + e~ annihilation interactions is a neutral off mass-shell photon or 
Z°, all events with only hadrons in the final state are the result of a pure quark- ant iquark pair 
creation process. Hard interactions in pp or pp collisions are much more complicated. Apart 
from the effect of the proton-remnant hadronisation, initial state QCD radiation predomi- 
nantly involving gluons is important. In addition gluon-gluon interactions dominate the hard 
scatter. Quark identification can only be performed statistically and in very limited kinematic 
areas. In e~p interactions at HERA there are fewer problems caused by the proton remnant 
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and, to first order, the virtual photon couples only to charged quarks. Electron kinematics 
alone can give the final state quark kinematics thus enabling a study of its hadronic fragments 
for direct comparison with quark fragmentation in e + e~ interactions. 

The momentum distribution of the hadrons from quark fragmentation in e + e~ interactions 
roughly scales with yfsZ, = E*, the overall centre of mass system (CMS) energy. The ratio of 
the momentum of a given charged hadron to the maximum energy that it could have, x p = 
2p hadrm / E* , is thus a natural variable in terms of which e + e~ experiments have described 
hadronic spectra. It has the advantage of being independent of the difficulties involved in any 
jet or thrust axis determination. The variable is manifestly not Lorentz invariant and it 
will therefore become necessary to identify an equivalent frame for e~p interactions in which 
comparisons can properly be made. 

The distribution D ± (x p ) = (1/N evts ) x dnf racks /dx p , being a charged track density nor- 
malised by the number of events, is termed the fragmentation function. It characterises the 
complete process which includes parton shower development as well as non-perturbative hadro- 
nisation. In principle the fragmentation function is unique only to a given quark and hadron 
species, but tagging techniques usually dictate that the data is effectively an integrated aver- 
age. It is a "soft" function rising rapidly as x p decreases but then turning over near x p = 0. 
A related distribution is that of the arbitrarily normalised logarithm of track energies where 
the turn-over becomes the "hump-backed" plateau (4|. 

The variable £ = ln(l/x p ) expands the turn-over region, and the expectation of the Modi- 
fied Leading Logarithmic Approximation (MLLA) to perturbative QCD predicts that for par- 
tons £ is distributed in a roughly Gaussian manner. Assuming Local Parton-Hadron Duality 
(LPHD) the same behaviour is expected for any type of hadron with the mean £ increasing as 
ln(E*) and the area under the Gaussian (multiplicity) increasing only slightly faster JJ. More 
detailed expectations of the fragmentation stem from the influence of quantum-mechanical 
coherence between gluons in the time-like parton shower |J. Effectively this reduces the 
available phase space for soft gluon emission to an angular ordered region. This coherence 
should then reduce the gradient of the energy dependence of the peak position. 

In e~p interactions there is no exact equivalent of the e + e~ centre of mass system. The 
hadronic (yp) CMS is one candidate. In the naive quark parton model (QPM) the final state 
then contains a quark which is back to back with the spectator system. This collaboration has 
already published |] xf = p z /Pmax distributions in this frame for the hemisphere opposite to 
the direction of travel of the incoming proton and outgoing remnant. 

A further boost along the (virtual photon)axis gives the Breit frame defined in this 
analysis such that the positive ^-direction is that of the incoming quark and the negative 
direction that of the incoming photon. This frame is defined such that the virtual photon 
is entirely space-like, having a momentum of p z = —Q. Within the QPM a quark has mo- 
mentum p z = +Q/2 before and p z = —Q/2 after scattering and therefore, in this frame, 
x p = 2p^ adron /Q is the equivalent of the e + e~ definition of a fragmentation variable^]. QCD 
radiation alters these values and implies that the collision is no longer collinear. In particular 

1 Note that since Q 2 is defined as a positive quantity one can refer to a value of 'Q' for the event. References 
discuss other possible variables in the Breit frame incorporating, for example, p» which gives directional 
information. Since the main aim of this paper is to compare with e + e _ experiments, we have chosen to remain 
with a compatible definition. 
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initial-state radiation of gluons has no equivalent in an e + e~ event. Note that whereas a given 
e + e _ experiment has one fixed value of E*, an e~p experiment results in a wide range of Q 
enabling evolution of the fragmentation function to be tested in a single experiment. Con- 
versely the fragmentation function should now be considered to be the result of a convolution 
with the structure functions of the individual quarks. This analysis will ignore this effect but 
will look for variations of fragmentation as a function of the event x as well as Q 2 . 

The Breit frame is defined solely by the kinematics of the virtual current and is concerned 
more with the quark-current rather than the proton-current interaction. The spectator system 
has reduced significance. Thus, despite some difficulties ||, the negative z hemisphere of this 
frame, rather than the 7p CMS, has been suggested as a better approximation to one half 
of the e + e~ CMS. Monte Carlo comparisons have been made to indicate the possibility 

of making the kind of measurements presented here. These analyses merely define a given 
hemisphere to contain current or target-associated tracks but in this analysis we shall also 
make some attempt to assess the confidence with which such associations may be made. 
Quantum-mechanically there is no absolute meaning for such associations, but as with the 
existence of jets, there are useful approximations which can be made. Since the jp CMS and 
the Breit frames are linked by a boost, they cannot in general both have the same negative z 
hemisphere content so the question of relative loss and contamination is worthy of empirical 
examination. 

2. The HI Experiment 



This analysis uses data from 1993 when the HERA e~p collider was run with a proton beam 
of 820 GeV and an electron beam of 26.7 GeV. A detailed description of the Hl-detector has 
been given elsewhere ||11|| . This paper will thus give only a short description of components 
vital for this analysis. 

Momentum measurements of charged tracks are provided, in the central region, by two 
cylindrical and co-axial drift chambers jn| for (r, <ft) measurement supplemented by two z- 



chambers; in the forward (proton) direction they are provided by three Radial and three 
Planar drift chamber modules |13| . These track detectors are inside a uniform 1.15 T magnetic 
field. Track segments from all devices are combined to give efficient detection and momentum 
measurement with 5p/p 2 better than 1% /GeV in the angular range used in this analysis, 
10° < 9 < 160°. A full Monte Carlo simulation based on GEANT [II gives an acceptable 
description of the effect of dead areas and reconstruction efficiency. 

In the polar angle range 4° < 9 < 153° the trackers are surrounded by a fine-grained liquid 
argon (LAr) sampling calorimeter |H| with lead and steel absorbers in the electromagnetic 



and hadronic sections, respectively. The calorimeter cells measure hadronic energy flow and 
the energy of the scattered electron for high Q 2 events. The LAr calorimeter is complemented 
by a backward electromagnetic, lead-scintillator, calorimeter (BEMC) covering the angular 
range 151° < 9 < 176°. The data are divided into two separate samples for low and high 
Q 2 where the scattered electron energy is well contained in either but not both of the two 
calorimeters. There are possible systematic energy scale errors taken to be ±1.7% for the 
BEMC and ±3% for the LAr calorimeter. 
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Data 


Non diffractive 


Good central 


Breit frame 


Sample 


events 


calorimetry 


energy flow 


12 < Q 2 < 80 GeV 2 


5927 


3141 


1945 


Q 2 > 100 GeV 2 


423 


373 


235 



Table 1: Size of data samples for this analysis. See text for definition of selections. 



3. Data Selection and Acceptance Corrections 



The source data for this analysis come from the late period of 1993 when magnetic field was 
reestablished in the HI detector. Selections to further limit data only to periods when there 
were extremely stable tracking conditions and when there was no coincident problem with 
electronic 'noise' in the HI calorimetry give a total integrated luminosity for this subset of 
« 150 nb _1 . 

A neutral current DIS event selection is made by demanding a well-reconstructed scattered 
electron of energy greater than 14 GeV. The selection 12 < Q 2 < 80 GeV 2 defines the BEMC 
sample, whereas the requirement for the scattered electron to have a polar angle in the range 
10° < 9 < 150° (resulting in the constraint Q 2 > 100 GeV 2 ) defines the LAr sample. These 
data sets are referred to as the low and high Q 2 samples, respectively. In terms of the usual 
DIS variables further cuts of W 2 > 3000 GeV 2 and, in the case of high Q 2 , y < 0.6 are 
applied to maintain high acceptance, and to avoid radiative corrections and contamination 
from photoproduction events. The final sample is better than 99.7% pure DIS. Additional 
cuts [0] are made to exclude possible diffractive events (about 6% of the sample) where the 
virtual photon may be said to interact with a Pomeron-like object [fl6 1 . There is no equivalent 



in the e + e~ interactions with which we intend to compare and such events are not adequately 
described by our DIS Monte Carlo generators with which we perform acceptance corrections. 
In practice this removal is done by demanding that there is at least 0.5 GeV deposit of 
observable energy in the polar region 4.4° < 9 < 15°. The resultant event numbers are given 
as the first column in Table [l|. 

Within these events, charged tracks are selected with Sp/p < 1 and with transverse mo- 
mentum p t > 150 MeV/c in the central tracker or momentum p > 500 MeV/c in the forward 
tracker. The tracks are also required to satisfy basic quality criteria on total number of 
digitisations, quality of vertex fits etc. Together these selections give a smooth variation in 
acceptance over the whole angular range. Cuts are made to remove tracks not originating 
from the primary interaction vertex but these fail to exclude roughly half of the ~ 4% of 
tracks arising from K® decays. 

Corrections for the effect of all of these cuts, and for the contribution of remaining K° 
and A decay products, are made by reconstructing Monte Carlo events with the same code 
and selections and then comparing these with results for primary charged particles at the 
generated event level. Commonly, LEPTO 6.1 [[[] is used with MRS H structure functions 
and with a Colour Dipole model |T/] for gluon radiation, but a version with matrix elements 



matched to parton showers in the MLLA approximation [Tj| has also been used and we note 



no significant change in any of the results. Since corrections are made bin-by-bin, in narrow 
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kinematic areas, and with ratios of generated to reconstructed track numbers the analysis is 
rather insensitive to Monte Carlo physics assumptions. 

4. Breit Frame Event Selection 

It is the aim of this analysis to select well-measured events requiring small acceptance correc- 
tions for analysis of tracks in the Breit frame of reference. This may be done using calorimetery 
information alone in order to avoid track reconstruction bias. First four-momentum vectors 
which correspond to the direction and energy of calorimeter clusters (assuming pion masses) 
are boosted to the Breit frame using the (x, Q 2 ) values calculated solely from the electron 
kinematics. The z < hemisphere, commonly referred to as the current fragmentation region, 
is then examined in more detail. To ensure that events have good and uniform track accep- 
tance a selection is made such that events are taken for further analysis if 95% or more of the 
observed energy in this hemisphere originates from the laboratory polar region 10° < 9 < 150°. 
This region is completely covered by the (LAr) calorimetry and thus has consistent calibra- 
tion. The numbers of such events surviving this cut are also given in Table Q in the column 
headed 'good central calorimetry'. 

For each event, the four-momenta of all energy clusters in the z < hemisphere are added 
and the resultant total energy E z<0 , is plotted as a fraction of the event Q against the resultant 
cosine of the Breit frame polar angle, cosQbf, in Fig. |T[ According to the QPM, a quark would 
be expected to have energy Q/2 and cosQ bf = — 1. At high Q 2 there is a tendency for the 
events to cluster in the area of this point. However, the spread of events is much greater at 
low Q 2 and there are a large number of events with small energy in the 'current' hemisphere 
and/or with a net cluster emerging at a large angle with respect to the beam axis. If, for 
a given event, the net energy flow is only just within the current hemisphere there could be 
some risk of bias in measuring the share of that energy. Simulated events samples behave in 
exactly the same way as the data and studies have shown that this spread is mainly because 
of reconstruction effects. Simple rejection of such events is a little problematic however as, at 
the generated level, some of these events show slightly different properties because of the loss 
of tracks to the target hemisphere through significant QCD radiation some of which would 
also have its equivalent in e + e~ interactions. We thus choose to test the possibilities of bias by 
performing the analysis both on all events and also by rejecting events below the line joining 
E z< o = 0, cosOfii? = —1 to E z<0 = Q/2, cosQbf = indicated in Fig. and referred to as 
the Breit frame energy flow selection. The event numbers surviving this topological selection 
are given as the final column in Table |l|. 

Having selected events using calorimetry, track properties in the Breit frame are now 
investigated, beginning with the distribution of cos8b, the cosine of the polar angle for charged 
particles, as shown in Fig. @(a,b). 

These data have bin-by-bin corrections for losses due to acceptance and inefficiencies cal- 
culated using simulated and reconstructed Monte Carlo events. For —1 < cos9b < 0.99, 
avoiding effects in this frame due to the beam pipe, the average correction factor is 1.35 for 
the low Q 2 sample and 1.42 for the high Q 2 sample. The correction factors are small and vary 
smoothly for cos# < 0.0. There is an obvious peak in the positive- z (beam) direction and also 
in the negative- z (quark) direction at high Q 2 but which is much less clear at low momentum 
transfer. 
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Figure 1: The total energy of the summed calorimeter cluster four momentum vectors in the 
z < hemisphere of the Breit frame is plotted as a fraction of the event Q against the cosine 
of the polar angle of the resultant vector, for (a) the low Q 2 and (b) the high Q 2 sample. The 
line indicates the cut for the secondary analysis described in the text. 




Figure 2: Distribution of the cosine of the Breit frame polar angle for tracks of (a) the low 
Q 2 and (b ) the high Q 2 sample, with statistical errors only. The open triangles show the data 
before the Breit frame energy flow selection. The solid line corresponds to the empirical fit 
described in the text, where the dashed line is the nominal quark contribution and the dotted 
line that of the target. 8 



In strict quantum-mechanical interpretation, no track can be said to 'belong' to either the 
jet or target remnant, wherever it lies in the Breit frame. It is clear from Fig. |^ however 
that the concept of current and target hemispheres is a good approximation as Q 2 increases 
but almost a question of mere definition at our lowest Q 2 values. We adopt two different 
procedures to estimate the possible loss or contamination involved in using a selection of 
one hemisphere of the Breit frame to represent the quark behaviour. In Fig. ^| there are 
empirical fits to Ae~~ a<Kl+x ^ a (1 — x) + Be~ b ^ 1 ~ x ' > ' 3 (1 + x), where x = cosOb- This and fits to other 
possible functional forms suggest that, in the Breit frame, the selection cos#b < 0.0 roughly 
equalises loss and contamination between the current and target fragmentation and that they 
are at the < 5% level at high Q 2 but at the ~ 20 — 25% level at low Q 2 . Our method of 
correction by Monte Carlo simulation accounts for migration between hemispheres caused by 
resolution effects. A separate study showed that these are dominated by measurement errors 
of the reconstructed electron giving an incorrect boost rather than track measurement errors. 
The magnitude of migration losses and contaminations are commensurate with the various 
empirical fits suggesting that these latter are indeed dominated by poor reconstruction. It is 
the approach of this analysis to adopt cos#b < 0.0 as our selection for tracks. If it is further 
required to associate these particle multiplicities with a quark we feel that an additional 
20(5)% systematic error on the particle multiplicity of low(high) Q 2 data should be imposed 
to describe our remaining ignorance of current /target separation. 

5. Fragmentation Functions 

The charged particle fragmentation functions, as defined in the introduction, are displayed in 
Fig. H(a,b) as a function of x p = 2p^ adron /Q for the current hemisphere of the Breit frame and 
for the two Q 2 intervals. These figures utilise the selections of the previous section including 
the Breit frame energy flow selection. The individual fragmentation functions, for positive 
and negative hadrons, are indistinguishable, at least with the present statistics in the current 
kinematic areas of investigation, and are treated inclusively for the rest of this paper. Changing 
the fragmentation variable to be £ = ln(l/x p ) and defining the fragmentation function to be 
D(£) = (1/N evts ) x dn trac ks/d^ results in the Gaussian shape of Fig. ffl(a,b) expected from 
the discussion in section 1 and seen also in e + e~ data. The area increases and the peak 
moves to higher values of £ at higher Q 2 . These distributions are also corrected bin- by-bin for 
acceptance and efficiency, although again the dependence is smooth and smallQ especially in 
the case of the variable £. Simulation studies show that the resolution in £ (5£ »s 0.12), again 
dominated by the uncertainties in the boost to the Breit frame is smaller than the bin-width 
of 0.2 in £ used throughout in this analysis. 

The low and high Q 2 data of Fig. [| have each been further subdivided to give a total 
of ten intervals of Q in order to study the evolution of the area, peak and width of the 
fragmentation function in more detail. The variation of the integrated area (equal to the 
mean charged particle multiplicity, < n ch >) is shown in Fig. |5| and the results are reproduced 
in Table 0. The systematic errors described at the end of section 4 dominate but are added in 
quadrature with statistical errors. The figure and table show the results obtained both with 
and without the Breit frame energy flow selections described in section 4 since there appears 

2 The effect of remaining QED radiative corrections to the Born term has been investigated using the 
Monte Carlo program Django 2.1 |l9[ . This analysis shows no effective variation over the peak area but an 
overall normalisation correction of +4.3 ± 0.3% has been included within our corrections, although it is small 
compared with other sources of error. 
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Figure 3: The fragmentation functions, D ± (x p ), for the current hemisphere of the Breit frame 
shown separately for positive and negative tracks, for (a) the low Q 2 and (b) the high Q 2 
sample with statistical errors only. 




Figure 4: The fragmentation functions for the current hemisphere of the Breit frame, D(£), 
for (a) the low Q 2 and (b) the high Q 2 sample, with statistical errors only and with simple 
Gaussian fits superimposed. 
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Figure 5: Average charged particle multiplicity in the current region of the Breit frame for 
data before (solid circles) and after (solid squares) the energy flow selection as a function of Q 
compared with a parameterisation (line) of one half of average corrected track multiplicities in 
e + e~ events as a function of E* . The results of a ZEUS analysis are shown as open squares, 
with points at the same Q being at different x. 

to be a significant difference, especially at low Q 2 . Since these data are corrected for known 
acceptance effects we ascribe this difference to QCD radiation depositing particles in the 
target hemisphere. It is seen similarly at the generated level of Monte Carlo. The difference is 
indicative of difficulties of detail inherent in a comparison with e + e~ data. A fit |2(J to many 
e + e~ data points^ is also shown for comparison as a function of E*. At low Q our results, 
and those of the ZEUS experiment [plj] , comparable to the HI results with no selection, show 
a significantly lower average charged multiplicity but compatible with differences that might 
be expected from differing QCD radiation effects. 

At a detailed level, agreement between particle multiplicities can also not be expected to 
be precise because of the comparative 'flavour democracy' expected in e + e~ annihilation. Jets 
initiated by 6-quarks, constituting some 22% of e + e~ events, have roughly 13% higher average 
charged multiplicity at LEP energies Such flavour effects are thus at the 3% level which 
is much less significant than our present errors and we make no attempt to correct for them. 
There is no doubt from this figure that Q is a suitable scaling quantity equivalent to E* in 
e + e~ interactions. Low energy, fixed target data indicated agreement between average 

3 Note that the published results refer to e + e~ events. Average track multiplicities have been reduced by a 
factor two to correspond to results for each quark and a further 8.1% (the average of available data) reduction 
has been made to account for K° and A decay tracks. 
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Q 2 Interval 


Energy Flow Selection 


Total Current Hemisphere 


(GeV 2 ) 


Multiplicity 


Peak 


Width 


Multiplicity 


Peak 


Width 


12 -> 15 


1.54±0.11 


1.41±0.13 


0.83±0.14 


1.35±0.08 


1.45±0.08 


0.72±0.08 


15 -> 20 


1.77±0.10 


1.47±0.10 


0.72±0.07 


1.63±0.08 


1.49±0.08 


0.82±0.09 


20 -> 40 


2.18±0.08 


1.63±0.09 


0.71±0.04 


2.02±0.06 


1.66±0.06 


0.74±0.04 


40 -> 60 


2.60±0.13 


1.91±0.11 


0.83±0.10 


2.39±0.11 


1.94±0.09 


0.91±0.11 


60 -> 80 


d.20±o.iy 


2.00±0.12 


0.83±0.12 


2.84±0.1o 


2.01±0.10 


0.87±0.13 


100 — > 175 


4.05±0.58 


2.23±0.39 


1.00±0.71 


3.71±0.50 


2.27±0.17 


0.75±0.27 


175 -> 250 


4.73±0.35 


2.12±0.34 


1.28±0.50 


4.05±0.29 


2.18±0.30 


1.26±0.48 


250 -> 450 


4.75±0.32 


2.72±0.29 


1.24±0.45 


4.00±0.25 


2.68±0.23 


1.09±0.32 


450 -> 1000 


5.43±0.39 


2.70±0.13 


0.80±0.15 


4.97±0.36 


2.71±0.16 


0.79±0.17 


1000 -> 8000 


6.27±1.20 


2.62±0.28 


0.88±0.39 


5.67±1.02 


2.63±0.26 


0.85±0.35 



Table 2: Average charged particle multiplicity, peak and width of the fragmentation function 
for the Q 2 intervals given using either the energy flow selection or the total current hemisphere 
of the Breit frame. 

charged particle multiplicities in the current hemisphere of the hadronic CMS frame and those 
of e + e~ with E* at the same W but later experiments Hl[p4f show significant disagreements 
at low Feynman-x. These analyses do not, in any case, attempt to demonstrate the relatively 
clear current-target separation seen, at least at high Q 2 , in this analysis. 

The peak and width (a) values of the fragmentation function may be found from Gaussian 
fits to the central area (± one unit of £ either side of the statistical mean) so as simultaneously 
to test MLLA predictions and to minimise dependence on Monte Carlo corrections which 
mainly affect the tails of the distributions. The simulated event studies referred to at the 
end of section 4 show that the resolution of a particle's polar angle in the Breit frame is 
much poorer than its momentum and indicate that at cos#b ~ 0, <5(cos#b) ~ ± 0.40(0.25) at 
low(high) Q 2 . This quantity does not enter into £ but affects peak and width measurements 
through uncertainty in the selection of the negative Breit frame hemisphere. As an estimator 
of the systematic error that might be introduced because of this selection we have repeated 
the above analyses using the cuts at the ± one sigma level of these resolutions and use extreme 
differences of results as an estimate of systematic errors. For the peak values, this estimate 
gives a result which is, in general, of the order of the statistical fitting errors. There is an 
additional source of systematic error of ~ ±6% in ^ pea k which arises from uncertainty in the 
boost due to the calorimeter energy scale errors given at the end of section 2. These errors 
dominate all other sources of systematic error coming, for example, from using different Monte 
Carlo generators for corrections. Width measurements are completely dominated by statistical 
errors. 

The results of these fits for both ^ pea k and ^ W idth variation as a function of Q are given in 
Table § and are shown in Fig. ^|(a,b), with all sources of error added in quadrature. The figure 
shows the results using our Breit frame energy flow selection but, as the results in the table 
indicate, there is only the suggestion of the peak values being slightly lower with the selection 
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rather than without. We may conclude that, within our errors, the parameters of this analysis 
are insensitive to the spectrum of hadronic fragments in QCD radiation lost to the target 
region and we therefore should still be able to compare with data from e + e~ interactions. The 
published e + e~ fragmentation function results have been fitted using the same procedure 
and are shown at the appropriate e + e~ E* values. Clear agreement is again evident using Q 
as the appropriate equivalent variable. The peak values are also in agreement with those of 
the ZEUS analysis [^TJ where no selection within the current hemisphere was made. 



If Q is taken to be equivalent to y/s^ and is normalised by an effective energy scale to give 
a variable Y = ln(Q/2A e jj) then, assuming gluon coherence, the predicted MLLA behaviour 
of the peak position and width is 

£ peak = 0.5Y + c 2 VY + O(l) 



iwidth = \jYij2cx 

where c\, c 2 are constants dependent only on the number of colours and flavours in QCD, 
Assets the scale of the mass of the final state fragmented hadrons, and 0(1) contains higher 
order corrections. Following || in assuming three flavours, we obtain A e jj= 0.24 ± 0.04 GeV 

and 0(1) = —0.38 ± 0.1 lQ in a combined fit (solid line) to the present peak and width data. 
These results are to be compared with the analysis @ of e + e~ i pea k evolution |f2]J^, ^5| , ^1 which 



gives A e jj= 0.21 ± 0.02 GeV and 0(1) = -0.32 ± 0.06. The dashed line of Fig. § shows the 
prediction at the generator level of the LEPTO 6.1 Monte Carlo with an assumption of gluon 



coherence in parton showers and using a Lund string model p7| for final hadronisation 



A straight line fit to ^ pea k as a function of ln(OJ gives a gradient of 0.75 ± 0.05. Reference 
claims a significant need for gluon coherence, noting that a gradient of unity would be the 
naive expectation if the multiplicity in a parton shower increased solely due to the constraints 
of longitudinal phase space. We find that abandoning the angular-ordering model of gluon 
coherence in the parton shower, but maintaining string hadronisation (dashed-dotted line in 
Fig. H|) gives almost indistinguishable evolution predictions. The results from HERWIG [28 



(not shown) with cluster hadronisation are similarly indistinguishable. A crude implementa- 
tion |29[ of an independent fragmentation model with or without (dotted line) gluon coherence 
gives effective gradients only slightly below unity (0.94 ±0.01 in this case, with our selections). 
However, we are aware of other implementations |30],|31| which impose energy conservation on 
hadronic distributions of transverse momentum that reproduce the same kind of 'string' ef- 
fects. We should conclude that the common e~p /e + e~ result does not constitute proof of the 
need for gluon coherence, and that the success of the MLLA-inspired model, even down to 



our lowest Q values, may be supporting claims El that the more fundamental cause of lower 



gradients is the increase of transverse momentum with energy. 

There is an expected effect on the peak position due to the limited opening angle of the jet 
implied by the selection of the current hemisphere in the Breit frame || . This systematically 
rejects lower x p , and hence higher £, tracks but there is a simultaneous inclusion of soft track 
contamination from the target. As we have shown, using one hemisphere of the Breit frame 

4 Thc evolution of the peak with Y is so nearly linear that there is a large correlation between the two 
parameters. This is resolved in our analysis by the width evolution which has no dependence on the O(l) 
term. 
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Figure 6: Evolution of (a) the peak position and b) the width (a) of a fitted Gaussian frag- 
mentation function compared with similar fits to e + e~ data and with the results of a ZEUS 
analysis, with points at the same Q being at different x. The solid line is a two-parameter 
simultaneous MLLA- expectation fit to the El data. See text for a discussion of the various 
model predictions. 
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Figure 7: Position of the fragmentation function peak as a function of Bjorken x, in the 
intervals (solid circles) 12 < Q 2 < 25 GeV 2 ; (open triangles) 25 < Q 2 < 45 GeV 2 ; (solid 
triangles) 45 < Q 2 < 80 GeV 2 ; and (open circles) 200 < Q 2 < 500 GeV 2 . The dashed lines 
refer to the ^ pea k expected from the fit at the relevant mean Q 2 value. 



roughly equalises these effects but the investigations |$3| of systematic errors which used cuts 
at the positions cosOb = ±0.40, ±0.25 as described earlier in this section show some evidence 
for a small systematic shift in the expected direction. The net effect is small and presently 
inconclusive. 

6. Dependence on x 

High Q 2 events are only accessible at high Bjorken x and this correlation is amplified by 
the W 2 and y cuts made in this analysis. Thus it is possible that the results presented so 
far are a reflected effect of a jet fragmentation evolution in terms of the parton momentum 
fraction x rather than Q 2 . The expectations of MLLA are that, although distributions in the 
target region should show some x dependence jL(J, the peak of £ in the current region should 



depend only on Q 2 . The possibility of such x dependence has been investigated, suffering 
the consequent loss of statistics, by defining a number of rectangular (x,Q 2 ) cells well within 
acceptance limits and somewhat broader in Q 2 than previously used. The corrected values of 
^ P eak are displayed in Fig. [7] as a function of x at constant Q 2 , with statistical errors only. The 
dashed lines refer to the expected values as interpolated from the fit to the function given in 
the previous section. It is clear that these data indicate no observable x dependence. 
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7. Conclusions 

This paper has shown that, when examined in the Breit frame as a function of Q, the basic 
features of the fragmentation function of charged particles in e~p interactions are close to 
those of quarks pair-created in e + e~ interactions when examined in their CMS as a function 
of yfsZ e . In short, to a good approximation quarks behave in the same way in both methods 
of production and Q is a suitable variable with which to describe the evolution of their 
fragmentation. On the contrary there is no evidence for any dependence on x. The data 
of this paper are consistent with the evidence from e + e~ data for the need to incorporate 
gluon coherence in the description of this fragmentation, but Monte Carlo models not making 
this assumption are also able to reproduce both sets of data quite well. 
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